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ABSTRACT 

We estimated the rate of comet and asteroid collisions with the terrestrial planets by calculating the orbits 
of 13000 Jupiter-crossing objects (JCOs) and 1300 resonant asteroids and computing the probabilities of 
. collisions based on random-phase approximations and the orbital elements sampled with a 500 yr step. 
[ — The Bulirsh-Stoer and a symplectic orbit integrator gave similar results for orbital evolution, but may give 



different collision probabilities with the Sun. A small fraction of former JCOs reached orbits with aphelia 
inside Jupiter's orbit, and some reached Apollo orbits with semi-major axes less than 2 AU, Aten orbits, 
, and inner-Earth orbits (with aphelia less than 0.983 AU) and remained there for millions of years. Though 

■ less than 0.1% of the total, these objects were responsible for most of the collision probability of former 

■ JCOs with Earth and Venus. We conclude that a significant fraction of near-Earth objects could be extinct 
comets that came from the trans-Neptunian region. 
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INTRODUCTION 

[ The main asteroid belt, the trans-Neptunian belt, and the Oort cloud are considered to be the main 
>^ ' sources of the objects that could collide with the Earth. Reviews of the asteroid and comet hazard were 
^ ; given by Ipatov (2000, 2001), and Bottke et al. (2002). Many scientists, e.g. Bottke et al. (2002), Binzel 
■ ■ ' et al. (2002), and Weissman et al. (2002), believe that asteroids are the main source of near-Earth objects 
(NEOs, i.e. objects with perihelion distance g<1.3 AU). Bottke et al. (2002) considered that there are 
200ibl40 km-sized Jupiter-family comets at q<1.3 AU, with ~80% of them being extinct comets. Duncan et 
al. (1995) and Kuchner (2002) investigated the migration of trans-Neptunian objects (TNOs) to Neptune's 
orbit, and Levison and Duncan (1997) studied the migration from Neptune's orbit to Jupiter's orbit. Ipatov 
and Hahn (1999) considered the migration of Jupiter-crossing objects (JCOs) with initial orbits close to the 
orbit of Comet P/1996 R2 and found that on average such objects spend about 5000 yr in orbits which cross 
both the orbits of Jupiter and Earth. Using these results and additional orbit integrations, and assuming 
that there are 5 x 10^ 1-km TNOs with 30<a<50 AU (Jewitt and Fernandez, 2001), Ipatov (2000, 2001) 
found that about 10^ 1-km former TNOs are Jupiter-crossers now and 10-20% or more 1-km Earth-crossers 
could have come from the Edgeworth-Kuiper belt into Jupiter-crossing orbits. In the present paper we use 
the estimates by Ipatov (2001), but now include a much larger number of JCOs. Preliminary results were 
presented by Ipatov (2002, 2003), who also discussed the formation of TNOs and asteroids. 

PROBABILITIES OF COLLISIONS OF NEAR-EARTH OBJECTS WITH PLANETS IN 
THE MODEL OF FIXED ORBITAL ELEMENTS 

As the actual collisions of migrating objects with terrestrial planets are rare, we use an approximation 
of random phases and orientations to estimate probabilities of collision for families of objects with similar 
orbital elements. We suppose that their semi-major axes a, eccentricities e and inclinations i are fixed, but 
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the orientations of the orbits can vary. When the orbit of a minor body crosses the orbit of a planet at a 
distance R from the Sun, the characteristic time to colhde, Tf, is a factor of k = v/vc = s/lajR — 1 times 
that computed with an approximation of constant velocity, where v is the velocity at the point where the 
orbit of the body crosses the orbit of the planet, and is the velocity for the same semi-major axis and a 
circular orbit. This coefficient k modifies the formulas obtained by Ipatov (1988a, 2000) for characteristic 
collision and close encounter times of two objects moving around the Sun in crossing orbits. These formulas 
depend also on the synodic period and improve on Opik's formulas when the semi-major axes of the objects 
are close to each other. As an example, at e=0.7 and a=3.06 AU, we have A;=2.26. 

Based on these formulas, wc calculated probabilities (1/T^) for ~1300 NEOs, including 343 Venus-crossers, 
756 Earth-crossers and 1197 Mars-crossers. The values of Tj (in Myr), /c, and the number of objects 
considered are presented in Table 1. We considered separately the Atens, Apollos, Amors, and several 
Jupiter- family comets (JFCs). The relatively small values of Tj for Atens and for all NEOs colliding with 
the Earth are due to several Atens with small inclinations discovered during the last three years. If we 
increase the inclination of the Aten object 2000 SG344 from i=0.1° to i=l°, then for collisions with the 
Earth we find T/=28 Myr and A;=0.84 for Atens and r/=97 Myr and A;=1.09 for NEOs. These times are 
much longer, and illustrate the importance of rare objects. Due to observational biases actual values of Tf 
can be greater than those in Table 1. 



Table 1. Characteristic collision times Tj (in Myr) of minor bodies with planets, coefficient 
k, and number Nj of simulated objects for the set of NEOs known in 2001. 
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ORBITAL EVOLUTION OF JUPITER-FAMILY COMETS AND RESONANT ASTEROIDS 

As the next step in estimating probabilities, we calculated the orbital evolution for thousands of test 
particles with initial orbits similar to known comets and asteroids, but having slightly different initial 
conditions. The results confirm that most of the collision probability comes from a handful of very rare 
cases in which the test particle is Earth-crossing for an extended period of time. 

For initial investigations of the migration of bodies under the gravitational influence of the planets, we 
used the integration package of Levison and Duncan (1994). In most cases we omitted the influence of 
Mercury and Pluto. Here and in Tables 2-3 and Figs. 1, 2 and 4 we present the results obtained by the 
Bulirsh-Stoer method (BULSTO code) with the integration step error less than eG[10~^-10~^], and in the 
next section we compare them with those of BULSTO at e<10^^^ and a symplectic method. 

In the first series of runs (denoted as nl) we calculated the evolution of 1900 JCOs moving in initial 
orbits close to those of 20 real JCOs with period 5<Pa<9 yr. In other series of runs, initial orbits were 
close to those of a single comet (2P, 9P, lOP, 22P, 28P, or 39P). For the 2P runs, we included Mercury in 
the integrations. We also investigated the orbital evolution of asteroids initially moving in the 3:1 and 5:2 
resonances with Jupiter. For the JCOs we varied only the initial mean anomaly ly. The number of objects in 
one run usually was <250. In most JCO cases the time r when perihelion was passed was varied with a step 
dT<l day (i.e., was varied with a step <0.2°). Near the r estimated from observations, wc used smaller 
steps. In most JCO cases the range of initial values of r was less than several tens of days. For asteroids, we 
varied initial values of v and the longitude of the ascending node from to 360° . The approximate values of 
initial orbital elements are presented in Table 2. We initially integrated the orbits for r5>10 Myr. After 10 
Myr we tested whether some of remaining objects could reach inside Jupiter's orbit; if so, the calculations 
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Table 2. Mean probability P=10~^Pr of a collision of an object with a planet (Venus=V, Earth=E, Mars=M) 
during its lifetime, mean time T (in Kyr) during which q<api, Tc=T/P (in Gyr), mean time Tj (in Kyr) spent 
in Jupiter-crossing orbits, mean time (in Kyr) spent in orbits with Q<4.2 AU, and ratio r of times spent in 
Apollo and Amor orbits. Results from BULSTO code at e~10~^-10~^. S denotes the sum for several series 
presented in the above lines. For iV=7349, 2P runs were excluded. 
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were usually continued. Therefore the results for orbits crossing or inside Jupiter's orbit were the same 
as if the integrations had been carried to Ufetimes of objects. For Comet 2P and resonant asteroids, we 
integrated until all objects were ejected into hyperbolic orbits or collided with the Sun. In some previous 
publications we have used smaller Ts, so these new data are more accurate. 

In our runs, planets were considered as material points so literal collisions did not occur. However, using 
the formulas of the previous section, and the orbital elements sampled with a 500 yr step, we calculated the 
mean probability P of collisions. We define P as Pt,/N, where Ps is the probability for all N objects of 
a collision of an object with a planet during its lifetime, the mean time T=Tj^/N during which perihelion 
distance q of an object was less than the semi-major axis api of the planet, and the mean time Tj during 
which an object moved in Jupiter-crossing orbits. The values of Pr=W^P, Tj and T are shown in Table 2. 
Here r is the ratio of the total time interval when orbits are of Apollo type (a>l AU, q=a{l — e)<1.017 AU) 
at e<0.999 to that of Amor type (1.017<g<1.3 AU) and Tc=T/P (in Gyr). In almost all runs T was equal 
to the mean time in orbits which cross the orbit of the planet and l/Tc was a probability of a collision per 
year (similar to 1/7/). The results showed that most of the probability of collisions of considered former 
JCOs with the terrestrial planets is due to a few objects each of which orbited for several Myr with aphelion 
(5<4.2 AU. Some had typical asteroidal and NEO orbits and reached Q<3 AU for several Myr. Large values 
of P for Mars in the nl runs were caused by a single object with a lifetime of 26 Myr. Total times spent 
by N JCOs and asteroids during their lifetimes in orbits typical for inner-Earth objects (lEOs, Q<0.983 
AU, Michel et al., 2000), Aten (a<l AU and Q>0.983 AU), A12 (g<1.017 AU and l<a<2 AU), Apollo, 
and Amor objects are presented in Table 3. These times for Earth-crossing objects were mainly due to a 
few tens of objects with high collision probabilities. Of the JCOs with initial orbits close to those of lOP 
and 2P, six and nine respectively moved into Apollo orbits with a<2 AU (A12 orbits) for at least 0.5 Myr 
each, and five of them remained in such orbits for more than 5 Myr each. The contribution of all the other 
objects to A12 orbits was smaller. Only one and two JCOs reached lEO and Aten orbits, respectively. 

One former JCO (Fig. la), which had an initial orbit close to that of lOP, moved in Aten orbits for 3.45 
Myr, and the probability of its collision with the Earth from such orbits was 0.344 (so Tc=10 Myr was even 
smaller than the values of Tf presented in Table 1; i.e., this object had smaller e and i than typical observed 
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Table 3. Times (in Myr) spent by N JCOs and asteroids during their lifetimes, with 
number of such objects in [ ]. Results from BULSTO code at e'-^10~^-10^*^. 
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Fig. 1. Time variations in a, e, q, Q, and i for a former JCO in initial orbit close to that of Comet lOP (a), 2P 
(b). For (a) at t<0.123 Myr Q>a>1.5 AU. Results from BULSTO code at e ~ lO"*^ - 10"^. 



Atens), greater than that for the 7850 other simulated former JCOs during their lifetimes (0.15). It also 

moved for about 10 Myr in lEO orbits before its collision with Venus, and during this time the probability 
Py=0.655 of its collision with Venus was greater (Py^S for the time interval presented in Fig. la) than 
that for the 7850 JCOs during their hfetimes (0.14). At t=0.12 Myr orbital elements of this object jumped 
considerably and the Tisserand parameter increased from J<3 to J>6, and J>10 during most of its lifetime. 
Another object (Fig. lb) moved in highly eccentric Aten orbits for 83 Myr, and its lifetime before collision 
with the Sun was 352 Myr. Its probability of collisions with Earth, Venus and Mars during its lifetime was 
0.172, 0.224, and 0.065, respectively. These two objects were not included in Table 2 except for the entry for 
A?^=7852. The data for Comet P/1996 R2 (line R2) were not included in the sums. Ipatov (1995) obtained 
the migration of JCOs into lEO and Aten orbits using the approximate method of spheres of action for 
taking into account the gravitational interactions of bodies with planets. The mean time Te during which a 
JCO was moving in Earth-crossing orbits is 9.6 x 10*^ yr for the 7852 simulated JCOs, and ~ 8 x 10^ yr for 
the 77,1 case. The ratio Ps of the number of objects colliding with the Sun to the total number of escaped 
(collided or ejected) objects was less than 0.015 for the considered runs (except for 2P). 

Some former JCOs spent a long time in the 3:1 resonance with Jupiter and with 2<a<2.6 AU. Other 
objects reached Mars-crossing orbits for long times. We conclude that JCOs can supply bodies to the regions 
which are considered by many scientists (Bottke et al., 2002) to belong to the main sources of NEOs. 
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Fig. 2. Distribution of migrating objects with their semi-major axes. The curves plotted in (b) at a=40 AU are 
(top-to-bottom) for sum, lOP, nl, 39P, 22P, 9P, 28P, and 2P. For Figs, (a) and (c), designations are the same. 
Results from BULSTO code at e ~ IQ-*^ - 10"^. 



In Fig. 2 we present the time in Myr during which objects had semi-major axes in the interval with a 
width of 0.005 AU (Figs. 2a-b) or 0.1 AU (Figs. 2c-d). At 3.3 AU (the 2:1 resonance with Jupiter) there is 
a gap for asteroids that migrated from the 5:2 resonance and for former JCOs (except 2P). 

For the nl data set, r7=0.12 Myr and, while moving in Jupiter-crossing orbits, objects had orbital periods 
P„<10, 10<Pa<20, 20<P„<50, 50<Pa<200 yr for 11%, 21%, 21%, and 17% of Tj, respectively. Therefore, 
there are three times as many JCOs as Jupiter-family comets (for which Pa<20 yr). Wc also found that 
some JCOs, after residing in orbits with aphelia deep inside Jupiter's orbit, transfer for tens of Myr to the 
trans-Neptunian region, either in low or high eccentricity orbits. We conclude that some of the main belt 
asteroids may reach typical TNO orbits, and then become scattered-disk objects having high eccentricities, 
and vice versa. The fraction of objects from the 5:2 resonance that collided with the Earth was only 1/6 of 
that for the 3:1 resonance. Only a small fraction of the asteroids from the 5:2 resonance reached a<2 AU. 

The distributions of migrating former JCOs and resonant asteroids in a and e and in a and i are presented 
in Fig. 4. For each picture we considered 250 migrating objects (288 for Fig. 4g-h), 100 intervals for a, 
and about 50-90 intervals for e and i. Different designations correspond to different number n of orbital 
elements (calculated with a step of 500 yr) which get in one bin (in Fig. 4 '<=' means <). All considered 
former JCOs very rarely reached low eccentricity orbits with 2<a<3.5 AU and ll<a<28 AU. There were 
many positions of objects when their perihelia were close to a semi-major axis of a giant planet, mainly of 
Jupiter. The pictures are different for different runs. 
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Fig. 3. Time variations in a, e, q, Q, and i for a former JCO in initial orbit close to that of Comet 2P. Results 
from a symplectic method at ds=10 days. 

COMPARISON OF ORBIT INTEGRATORS 

To determine the effect of the choice of orbit integrators and convergence criteria, we made additional runs 
with BULSTO at £=10"^^ and £=10-^^ and with a symplectic integrator. The orbital evolution of 5400 
JCOs was computed with the RMVS3 code. For the symplectic method we used an integration step dg of 3, 
10, and 30 days. We find that for the purposes of this paper, the differences between integrator choices (at 
ds<10 days) are comparable to the differences between runs with slightly different initial conditions. Our 
interpretation is that 1) very small numbers of particles contribute most of the collision probabilities with 
the terrestrial planets, 2) runs with larger numbers of particles are more reliable, and 3) small differences 
in initial conditions or in the errors of the orbit integrators modify the trajectories substantially. 

To illustrate these points, Tables 4-5 present the results obtained by BULSTO at s<10~^^ and the 
symplectic method at ds<lO days. Most of the results obtained with these values of s and dg are statistically 
similar to those obtained for 10-^<e<10-^. For example, a few objects spent millions of years in Earth- 
crossing orbits inside Jupiter's orbit (Figs. 1 and 3), and their probabilities of collisions with the Earth 
were thousands of times greater than for more typical objects. For series nl with RMVS3, the probability 
of a collision with Earth for one object with initial orbit close to that of Comet 44P was 88.3% of the total 
probability for 1200 objects from this series, and the total probability for 1198 objects was only 4%. This 
object and the object presented in Fig. 3 and in the first line of Table 4 were not included in Table 5 with 
A^=1199 for nl and with A^=250 for 2P, respectively. For the resonance 3:1 at ^5=10 days, 142 objects spent 
140 and 84.5 Myr in lEO and Aten orbits, respectively, even longer than for £~10-^-10-^. Additionally, up 
to 40 Myr and 20 Myr were spent in such orbits by two other objects which had estimated probabilities of 
collisions with the terrestrial planets greater than 1 for the calculated sets of orbital elements. For the 2P 
runs at e<10~^^ and A'^=100, the calculated objects spent 5.4 Myr in Apollo orbits with a<2 AU. 



Table 4. Times (Myr) spent by three objects in various orbits, and probabilities of collisions with 
Venus {py), Earth (pe), and Mars {pm) during their lifetimes Tu (in Myr). Symplectic method. 
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The values of Pr presented in Table 5 are usually of the same order of magnitude as those in Table 2, and 
the difference between the data presented in these tables is comparable to the differences between different 
runs belonging to a series. For Earth and Venus, the values of Pr presented in both tables are about 1-4 for 
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ig. 4. Distribution of migrating objects in semi-major axes, eccentricities, and inclinations for objects in initial 
rbits close to that of lOP (a-b,e), 2P (c-d.f), and the resonance 3:1 with Jupiter (g-h). BULSTO code at 
~ 10~^ — 10"*^. For (e-f,h) it was considered that an object disappeared when perihelion distance became less 
ban 2 radii of the Sun. In other cases, objects disappered when they collided with the Sun. 
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9P, 22P, 28P and 39P. For 28P and 39P with the symplectic method Pr is about twice that for BULSTO. 
For lOP, Pr is several times larger than for the above series, and for 2P it is several times larger than for lOP. 
For the nl run, Pr>4 for Earth. The ratio of Pr to the mass of the planet was typically several times larger 
for Mars than for Earth and Venus. The main difference in Pr was found for the 3:1 resonance. In this case 
greater values of Pr were obtained for ds=W days. As noted above, a few exceptional objects dominated 
the probabilities, and for the 3:1 resonance two objects, which had collision probabilities for calculated sets 
of orbital elements greater than unity for the terrestrial planets, were not included in Table 5. These two 
objects can increase the total value of Pr for Earth by a factor of several. 

In Tables 2 and 5 wc present the mean time (in Kyr) spent in orbits with Q<4.2 AU. It can differ by 
three orders of magnitude for different series of runs. In Table 5 in [ ] we present the number of objects 
each of which got Q<4.2 AU during at least 1000 yr, the second number AT^ in [ ] means the same but only 
for g<1.017 AU. For series 2P and the 3:1 resonance with Jupiter both values of were almost the same, 
for lOP they differed by a factor of 2.9 (symplectic runs) or 3.7 (BULSTO), and for other series only a small 
portion of decoupling objects crossed the orbit of the Earth. For most runs (except for 2P and asteroids) 
the number of objects which got Q<4.7 AU was several times larger than that for Q<4.2 AU. 

We also did some symplectic runs with ds=30 day. For most of the objects we got similar results, but 
about 0.1% of the objects reached Earth-crossing orbits with a<2 AU for several tens of Myr and even lEO 
orbits. These few bodies increased the mean value of P by a factor of more than 10. With ^^=30 days, four 
objects from the runs nl, 9P, lOP had a probability of collisions with the terrestrial planets for calculated 
sets of orbital elements greater than 1 for each, and for 2P there were 21 such objects among 251 considered. 
For resonant asteroids, we also obtained much larger values than those with BULSTO for P and T with 
RMVS3 at ds=SO days, and similarly for the 3:1 resonance even at ds=10 days. For this resonance it may 
be better to use ds<10 days. Probably, the results of symplectic runs at ds=SO days can be considered as 
such migration that includes some nongravitational forces. 



Table 5. Probabilities of collisions with the terrestrial planets. Designations are same as those for Table 2. 
is the mean time (in Kyr) spent in orbits with Q<A.2 AU; in [ ] we present the number TV,; of objects each of 
which got Q<A.2 AU during at least 1000 yr, the second number in [ ] means the same but only for g<1.017 
AU. Results from the BULSTO code at e<l{)~^^ and a symplectic method at (jg<10 days. 
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313 


3.32 


10.0 


324 


32.3 


35. 


0.29 


585(250/250] 


9P 


10'^ 


400 


1.37 


3.46 


2.53 


3.26 


7.84 


2.40 


1.62 


23.8 


14.7 


1.1 


128 


8.0(13/3] 


lOP 


<10'^ 


450 


14.9 


30.4 


2.04 


22.4 


41.3 


1.84 


6.42 


113. 


17.6 


1.5 


85 


44.(70/24] 


22P 


lO'^ 


250 


0.68 


2.87 


4.23 


1.39 


4.96 


3.57 


0.60 


11.5 


19.2 


1.5 


121 


0.6(3/0] 


28P 


10*^ 


250 


3.87 


35.3 


9.12 


3.99 


59.0 


14.8 


0.71 


109. 


154 


2.2 


535 


3.3(7/0] 


39P 


10<^ 


250 


2.30 


2.68 


1.17 


2.50 


4.22 


1.69 


0.45 


7.34 


16.3 


2.2 


92 


0.5(2/0] 


3 : 1 


< 10-12 


70 


1162 


1943 


1.67 


1511 


5901 


3.91 


587 


803 


1.37 


4.6 


326 


8400(70/70] 


3 : 1 


lO'^ 


142 


27700 


8617 


0.31 


2725 


9177 


3.37 


1136 


9939 


8.75 


16. 


1244 


5000(142/140] 


5 : 2 


10-12 


50 


130 


113 


0.87 


168 


230 


1.37 


46.2 


507 


11.0 


1.4 


166 


512(50/4] 


5 : 2 


10"^ 


144 


58.6 


86.8 


1.48 


86.7 


174 


2.01 


17. 


355 


20.9 


1.7 


224 


828(144/13] 



In the case of close encounters with the Sun (Comet 2P and resonant asteroids), the probability Ps of 
collisions with the Sun was larger for RMVS3 than for BULSTO, and for 10-i^<£:<10-i2 it was greater 
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than for 10 ^<£<10 ^ (P5=0.75 for the 3:1 resonance at ds=3 days). This probabihty is presented below 
for several runs: 



Table 6. Probability of collisions with the Sun 





e = 10-^^ 


e = 10-''' 


e = 10-y 


£ = 10-« 


ds = 10 days 


ds = 30 days 


Comet 2P 


0.88 


0.88 


0.38 


0.32 


0.99 


0.8 


resonance 3 : 1 


0.46 


0.5 


0.156 


0.112 


0.741 


0.50 


resonance 5 : 2 




0.06 


0.062 


0.028 


0.099 


0.155 



For Comet 2P the values of Tj were much smaller for RMVS3 than those for BULSTO and they were 
smaller for smaller e; for other runs these values do not depend much on the method. In the direct 
integrations reported by Valsecchi et al. (1995), 13 of the 21 objects fell into the Sun, so their value of 
^5=0.62 is in accordance with our results obtained by BULSTO; it is less than that for e=10~^^, but 
greater than for e=10~^. Note that even for different Ps the data presented in Tables 2 and 5 usually are 
similar. As we did not calculate collision probabilities of objects with planets by direct integrations, but 
instead calculated them with the random phase approximation from the orbital elements, we need not make 
integrations with extremely high accuracy. Ipatov (1988b) showed that for BULSTO the integrals of motion 
were conserved better and the plots of orbital elements for closely separated values of e were closer to one 
another at 10~^<e<10~^. The smaller the value of e, the more integrations steps are required, so £:<10~^^ 
for large time intervals are not necessarily better than those for 10~^<£:<10~^. Small e is clearly necessary 
for close encounters. Ipatov and Hahn (1999) and Ipatov (2000) found that former JCOs reached resonances 
more often for BULSTO than for RMVS3 at ds=SO days. Therefore we made most of our BULSTO runs 
with 10~^<e<10~^. For a sjnnplectic method it is better to use smaller dg at a smaller distance R from the 
Sun, but in some runs R can vary considerably during the evolution. 

W. Bottkc pointed out that H. Lcvison showed that it is difficult to detect solar collisions in any numerical 
integrator, so he removed objects with qKqmin- The results presented above were obtained considering 
collisions with the Sun, but we also investigated what happens if we consider gmin equal to ks radii rs of the 
Sun. For ks=2, some results are presented in Fig. 4e-f,h. The only difference with the runs that considered 
collisions with the Sun is that for those runs for series 2P and lOP and for the 3:1 resonance, some objects 
reached 90°<z<180° (mainly with 2<a<3.5 AU) (Fig. 4b,d,g). For ks=2 there were no comets with i>90° 
and there were only a few orbits of asteroids with i>90° (Fig. 4e-f,h). The consideration of gmin at ks=3 
did not influence the collision probabilities with the terrestrial planets or getting orbits with a<2 AU. For 
example, with BULSTO for the two objects with the largest collision probabilities, the time spent in orbits 
with a<2 AU decreased by only 0.3% for 2P at ^5=3 and was the same for lOP at ^5=10. 

MIGRATION FROM BEYOND JUPITER TO THE TERRESTRIAL PLANETS 

According to Duncan et al. (1995), the fraction Ptnj of TNOs reaching Jupiter's orbit under the influence 
of the giant planets in 1 Gyr is 0.8-1.7%. As the mutual gravitational influence of TNOs can play a larger 
role in variations of their orbital elements than collisions (Ipatov, 2001), we considered the upper value of 
Ptnj- Using the total of 5 x 10^ 1-km TNOs with 30<a<50 AU, and assuming that the mean time for 
a body to move in a Jupiter-crossing orbit is 0.12 Myr, we find that about A^jo=10'^ 1-km former TNOs 
are now Jupiter-crossers, and 3000 are JFCs. Using the total times spent by N simulated JCOs in various 
orbits, we obtain the following numbers of 1-km former TNOs now moving in several types of orbits: 





method 


series 


lEOs 


Aten 


A12 


Apollo 


Amor 


3100 


BULST0+RMVS3 


nl 








480 


1250 


900 


1900 


BULSTO 


nl 








25 


720 


1000 


7350 


BULSTO 


all without 2P 


120 


40 


250 


2500 


1750 


7852 


BULSTO 


all 


110 


950 


4500 


7200 


1880 
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For example, the number of lEOs NiEOs=^JotiEo/i^jij), where tjEO is the total time during which Nj 
former JCOs moved in lEOs' orbits, and Njtj is the total time during which Nj JCOs moved in Jupiter- 
crossing orbits. As we considered mainly the runs with relatively high migration to the Earth, the actual 
number of NEOs is smaller by a factor of several, the actual portion of lEOs and Atens can be smaller 
and that for Amors can be larger than those in the lines in the Table at iV>7000. Even if the number of 
Apollo objects is smaller than that based on nl runs, it may still be comparable to the real number (750) 
of 1-km Earth-crossing objects (half of them are in orbits with a<2 AU), although the latter number does 
not include those in highly eccentric orbits. 

The values of the characteristic time (usually Tc) for the collision of a former JCO or a resonant asteroid 
with a planet (see Tables 2 and 5) are greater than the values of Tf for NEOs in Table 1, Tc~ 1.1 Gyr for 
7852 objects, so we expect that the mean inclinations and eccentricities of unobserved NEOs are greater 
than those for the NEOs that are already known. Jedicke et al. (2003) found similar results. On average, 
the values of Tc for our nl series and for most of our simulated JCOs were not greater than those for our 
calculated asteroids, and migrating Earth-crossing objects had similar e and i for both former JCOs and 
resonant asteroids. Former JCOs, which move in Earth-crossing orbits for more than 1 Myr, while moving in 
such orbits, usually had larger P and smaller e and i (sometimes similar to those of the observed NEOs, see 
Figs. 1 and 3) than other JCOs. It is easier to observe orbits with smaller values of e and i, and probably, 
many of the NEOs moving in orbits with large values of e and i have not yet been discovered. About 1% 
of the observed Apollos cross Jupiter's orbit, and an additional 1% of ApoUos have aphelia between 4.7-4.8 
AU, but these Jupiter-crossers arc far from the Earth most of time, so their actual fraction of ECOs is 
greater than for observed ECOs. The fraction of Eartli-crossers among observed Jupiter-family comets is 
about 10%. This is a little more than T/Tj for our nl runs, but less than for 7850 JCOs. For our former 
resonant asteroids, Tj is relatively large (RiO.2 Myr), and such asteroids can reach cometary orbits. 

Comets are estimated to be active for Taef^lO^-lO^ yr. Tact is smaller for closer encounters with the Sun 
(Weissman et al., 2002), so for Comet 2P it is smaller than for other JFCs. Some former comets can move 
for tens or even hundreds of Myr in NEO orbits, so the number of extinct comets can exceed the number of 
active comets by several orders of magnitude. The mean time spent by Encke-type objects in Earth-crossing 
orbits is >0.4 Myr (even for c/min)- This time corresponds to >40-400 extinct comets of this type. Note that 
the diameter of Comet 2P is about 5-10 km, so the number of smaller extinct comets can be much larger. 

The above estimates of the number of NEOs are approximate. For example, it is possible that the number 
of 1-km TNOs is several times smaller than 5 x 10^, while some scientists estimated that this number can 
be up to 10^^ (Jcwitt, 1999). Also, the fraction of TNOs that have migrated towards the Earth might be 
smaller. On the other hand, the above number of TNOs was estimated for a<50 AU, and TNOs from more 
distant regions can also migrate inward. Probably, the Oort cloud could also supply Jupiter-family comets. 
According to Asher et al. (2001), the rate of a cometary object decoupling from the Jupiter vicinity and 
transferring to an NEO-like orbit is increased by a factor of 4 or 5 due to nongravitational effects (see also 
Fernandez and Gallardo, 2002). This would result in larger values of Pr and T than those shown in Tables 
2 and 5. Our estimates show that, in principle, the trans-Neptunian belt can provide a significant portion 
of the Earth-crossing objects, although many NEOs clearly came from the main asteroid belt. It may be 
possible to explore former TNOs near the Earth's orbit without sending spacecraft to the trans-Neptunian 
region. According to our results, many former Jupiter-family comets can have orbits typical of asteroids, 
and collide with the Earth from typical NEO orbits. 

Based on the collision probability P = 4 x 10~^ we find that 1-km former TNOs now collide with the 
Earth once in 3 Myr. This value of P is smaller than that for our nl runs, does not include the 'champions' 
in collision probability, and is only 1 /20 of that for our 7852 JCOs. Using P = 4 x 10"^ and assuming the 
total mass of planetesimals that ever crossed Jupiter's orbit is ^lOOrn®, where is the mass of the Earth 
(Ipatov, 1993, 2000), we conclude that the total mass of bodies that impacted the Earth is 4 x 10~^m®. If 
ices comprised only half of this mass, then the total mass of ices Mice that were delivered to the Earth from 
the feeding zone of the giant planets is about the mass of the terrestrial oceans (~ 2 x lO^^me). 

The calculated probabilities of collisions of objects with planets show that the fraction of the mass of the 
planet delivered by short-period comets can be greater for Mars and Venus than for the Earth (compare 
the values of P/rUpi using P from Tables 2 and 5, where rupi is the mass of the planet). This larger mass 
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fraction would result in relatively large ancient oceans on Mars and Venus. On the other hand, there is the 
deuterium/hydrogen paradox of Earth's oceans, as the D/H ratio is different for oceans and comets. Pavlov 
et al. (1999) suggested that solar wind-implanted hydrogen on interplanetary dust particles could provide 
the necessary low-D/H component of Earth's water inventory. 

Our estimate of the migration of water to the early Earth is in accordance with Chyba (1989), but is 
greater than those of Morbidelli et al. (2000) and Levison et al. (2001). The latter obtained smaller values 
of Mice, and we suspect that this is because they did not take into account the migration of bodies into orbits 
with Q<4:.2 AU and q<l AU. Perhaps this was because they modeled a relatively small number of objects, 
and Levison ct al. (2001) did not take into account the influence of the terrestrial planets. In our runs the 
probability of a collision of a single object with a terrestrial planet could be much greater than the total 
probability of thousands of other objects, so the statistics are dominated by rare occurrences that might 
not appear in smaller simulations. The mean probabilities of collisions can differ by orders of magnitude for 
different JCOs. Other scientists considered other initial objects and smaller numbers of JCOs, and did not 
find decoupling from Jupiter, which is a rare event. We believe there is no contradiction between our present 
results and the smaller migration of former JCOs to the near-Earth space that was obtained in earlier work, 
including our own papers (e.g. Ipatov and Hahn, 1999), where we used the same integration package. 

From measured albedos, Fernandez et al. (2001) concluded that the fraction of extinct comets among 
NEOs and unusual asteroids is significant (at least 9% are candidates). Rickman et al. (2001) believed 
that comets played an important and perhaps even dominant role among all km-size Earth impactors. In 
their opinion, dark spectral classes that might include the ex-comets are severely undcrreprcscntcd (sec also 
Jewitt and Fernandez, 2001). Our runs showed that if one observes former comets in NEO orbits, then it 
is probable that they have already moved in such orbits for millions (or at least hundreds of thousands) 
years, and only a few of them have been in such orbits for short times (a few thousand years). Some former 
comets that have moved in typical NEO orbits for millions or even hundreds of millions of years, and might 
have had multiple close encounters with the Sun, could have lost their mantles, which causes their low 
albedo, and so change their albedo (for most observed NEOs, the albedo is greater than that for comets; 
Fernandez et al., 2001) and would look like typical asteroids. Typical comets have larger rotation periods 
than typical NEOs (Binzel et al., 1992), but, while losing considerable portions of their masses, extinct 
comets can decrease these periods. In future we plan to consider a larger initial number of objects initially 
located beyond Jupiter in order to better estimate their probabilities of migration to a near-Earth space. 

CONCLUSIONS 

Collision statistics for the terrestrial planets are dominated by very small numbers of bodies that reach 
orbits with high collision probabilities, so it is essential to consider very large numbers of particles. The 
initial conditions for the orbit integrations appear to matter more than the choice of orbit integrator. Some 
Jupiter-family comets can reach typical NEO orbits and remain there for millions of years. While the 
probability of such events is small (about 0.1% in our runs and perhaps smaller for other initial data), 
nevertheless the majority of collisions of former JCOs with the terrestrial planets are due to such objects. 
The amount of water delivered to the Earth during planet formation could be about the mass of the Earth 
oceans. From the dynamical point of view there could be (not 'must be') many extinct comets among the 
NEOs. For better estimates of the portion of extinct comets among NEOs we will need orbit integrations 
for many more TNOs and JCOs, and wider analysis of observations and craters. 
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